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Abstract 

The present paper examines the active conttol of parametric resonance in 
modified Rayleigh-Duffing oscillator. We used the method of averaging to 
obtain steady-state solutions. We have found the critical value of the para- 
metrical amplitude which indicates the boundary layer where the control is 
efficient in reducing the amplitude vibration. We find also the effects of 
excitation parameters and time-delay on dynamical of this system with the 
principal parametric resonance. We obtain also for this oscillators the Hopf 
bifurcation or saddle-node bifurcation for certains values of parametric pa- 
rameters and time-delay and we have studied the influence of parameter k2 
which is one parameter which modify the ordinary Rayleigh-Duffing oscilla- 
tor. We have discussed the appropriate choice of the time-delay and control 
gain. We finally studied the stability of fixed point and it is found that the 
appropriate choice of the time-delay can broaden the stable region of the 
non-trivial steady-state solutions enhance the control performance. Numeri- 
cal simulations are performed in order to confirm the analytically results. 

keywords: Active control, parametric resonance, primary, modified Rayleigh- 
Duffing oscillator, stability, Hopf burfication, amplitude vibration. 



1 Introduction 

In recent years, a twofold interest has attracted theoretical, numerical, and exper- 
imental investigations to understand the behavior of nonlinear oscillators. The- 
oretical (fundamental) investgation reveals their rich and complex behavior, and 
the experimental (self-excited oscillators) describes the evolution of many biologi- 
cal, chemical, physical, mechanical, and industrial systems |[Tl|2j[3l. Recently, the 
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chaotic behavior of these oscillators is exploited in the field of communication for 
coding information [4]. The considerable efforts have been devoted to the control 
of vibrating structures in various fields of fundamental and applied sciences. Two 
major aims in the scope of the researchers are: reduce the amplitude of vibrations, 
inhibit chaos and escape from potential well. Applications are well encountered 
in structural mechanics (see Refs. [1-5]). Parametric excitation occurs in a wide 
variety of engineering application (Refs. [10-14]). Among various types of control 
strategies, the active feedback control usually developed by means of electromag- 
netic force or by a mechanical device ( Ref.[5-9]) and the parametric. In liTSll 
Xin-Ye and al demonstrated that in the dynamical behaviour of a parametrically 
excited Duffing- Van der Pol oscillator under linear-plus-nonlinear state feedback 
control with a time-delay, nontrivial steady state responses may lose their stability 
by a saddle-node or Hopf bifurcation as parameters vary. M.Siewe Siewe and al 
|[T6l studied the dynamics of a parametrically excited Rayleigh-Duffing oscillator 
state feedback control with a time-delay, nontrivial steady state solutions may lose 
their stability by a saddle-node or Hopf bifurcation as parameters vary. 

Our objective of this paper is to control the amplitude vibration of a paramet- 
rically excited modified Rayleigh-Duffing oscillator with time-delayed feedback 
position and linear velocity terms. The modified Rayleigh-Duffing oscillator that 
we consider in this work is governed by following equation: 

X -\- eax + e(3x'^ + ejx^ — 2ekixx + ek2X^x + [5 + ea cos 2ujot)x 
+eXx^ = 0. (1) 

We used the method of averaging for study this oscillator. The amplitude peak of 
the parametric resonance can be reduced by means of a correct choice of the time- 
delay and the feedback gain. We discussed how the existence region of steady-state 
solutions is modified by the feedback control and we showed the existence regions 
of the nontrivial solution in the plane of the parametric excitation amplitude and the 
detuning parameter for an uncontrolled system, a controlled system without time- 
delay, and those with time-delays corresponding to the minimum and maximum 
values of an appropriate equivalent damping. A bifurcation analysis and parametric 
excitation-response and frequency-response curves are presented. We found the 
effect of k2 on the control. Then we performed a stability analysis of the bifurcation 
of the model and we derived sufficient conditions for stable nontrivial solutions in 
order to exclude the presence of modulated motion. 

2 Model and linear control 
2.1 Model 

Consider the modified Rayleigh-Duffing oscillator equation which describes the 
glycolytic reaction catalized by phosphofructokinase, namely the Selkov equations 
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and abstract trimolecular chemical reaction ifTTl . 

Perturbing this system by Duffing force and parametric excitation force, we 
obtained the parametric dissipative modified Rayleigh-Duffing oscillator which is 
expressed by eq. ([T]). We consider a single-degree-of-freedom model with a nonlin- 
ear soft spring, a linear damper and a linear displacement feedback control system. 
The model is described by the following differential equation system: 



X + eax + ef5x'^ + eji"^ — 2ekixx + ek2X^x + {6 + ea cos 2ujQt)x 
+€Xx^ - ehx{t - r) - ecx{t - r) = 0, (3) 

where x,x, x are displacement, velocity and acceleration respectively; ujq is the 
natural frequency, e is a positive parameter, a, /?, 7, ki,k2, 6, A, a are constants , 
b and c are the control gain parameters and r is time-delay parameter. In this 
equation, —ebx{t — r) — ecx{t — r) is the feedback force in the system. The 
resonances appear when 6 = uiq + e6i. 

2.2 Linear control 

In this part, we used the method of averaging to study the eq. Q and to research 
the controlling domain. In the absence of the parameter e, eq.Q reduces to 

X + WqX = 0, (4) 

with the solution: 

X = A cos{ujot + (j)). (5) 

If e ^ 0, A and </> depend on time t. In this condition, differentiating ^ and 
using the method of averaging, we found that the eq. Q is divided into: 

dA 

— — = —eFlA cosfwoi + 4')^ sinfwot + sinfwoi + ^fy (6) 
at 



A— = —eFlA cosfwoi + 4')^ ~^ sinfwot + 0)1 cosfcjoi + 1;^), (V) 
at 



where 



F = F{x, x) = —ax — fiiP' — 73;^ + 2kixx — k2x'^x 

— {5i + a cos 2u}ot)x — Xx^ + bx{t — r) + cx{t — r). (8) 

When applying the hypothesis of the method of averaging (^4 and (p are small 
for one period for small e), we obtain 

rTo 



dA 




Itt 


To I 


dS 


-rf 


dt 


Jo Jo 



A- 
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cos{uJot + (j))dt, 



(10) 



where Tq is the natural period of this oscillator. Using x and his derivative, eq.([8]) 
becomes 

F[Acos{ijjQt + (j)), —Asm{uJot + cp)] = aujoAsm{u)Qt + (/>) 

-/3ojqA'^ sin^iojot + </>) + 7^o^^ sin^(a;ot + 0) + 2kiUJoA'^ sm{ojot + (f>)x 

cos{ujot + (j)) — k2A^u>Q sin^ (ujQt + (j)) cos(a;ot + </>) 

— {6i + acos2a;ot)Acos(wot + </>) — Aj4^ cos^(a;oi + 0) 

+bAcos[uJo{t — t) + (p] — cujoAsm[uJo{t — r) + i;^], (11) 

where A = A{t — r) and (j) = (l){t — r). 

Substituting ( [TT] ) into the eq.(|9]) and eq.dOl and evaluating the integrals, we 

find 



dA 
~dt 



-^[uJoA{2a + ^jujqA^) - aAsin'. 



+2bA sin( 



+ loqt) — 2cu>oAcos{<p — (p + ojot)], (12) 



A- 



dt 



I [25iA + ^{k2 + 3X)A^uj^ + aA cos 2</> 



cosf 



+ Wo''") + 2cujoA sin( 



+ a;or)]. (13) 



Eqs.( 12 1 and ( 13 1 show that A and (p are 0(e). Expanding in Taylor series A 

A = A{t-T) = A{t)-TA{t)+T'^A{t) + 

<p = (p{t-T) = (P{t)-T^{t)+T^4>{t) + .... 



(14) 
(15) 



Eqs.(14l and (15 1 indicate that we can replace A and (p by A and (p in eq.(12) 
and (13 1 since A, (/> and A,cp are 0(e) and O(e^) respectively. This reduces an 



infinite-dimensional problem in functional analysis to a finite-dimensional problem 
by assuming that the product er is small lfT6l . By setting t = i/e as the new time 
scale, we have the following averaged equations: 



A 



A(P 



A, 



[ujQ{—2a — -^ujqA ) + a sin 20 — 26sincjo''" 



4 L vv 2 
+2cujQ cos ojqt], 
A 1 

— [26i + - (/c2 + 3A)^^Wo + « cos 2(1) 
+2ccjo sin ojqt]. 



26 cos ujqt 



(16) 



(17) 



Seeking steady state, we have ^4 = and Acp = in ( 16 1 and ( 17 1. A trivial 



solution is A = and we find that the non trivial solutions verify the equation: 



Al - 2PAl + Q 



0, 



(18) 
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where the quantities P and Q are given by 



1 3 o 

-[-^0OQ{—2u}Qa — IbsinujQT + 2cujq cosujqt) 



|7M + 3(fc2 + 3A)2a;4^2 
2^0(^2 + 3A)(2(5i — 26cosa;oT + 2cuJo sinwo''")]) (19) 



Q = Q ^ « -—^[{-2ujoa-2bsmujoT + 2cuJocosuJoTf+ 

iT^^o + 4(^2 + 3A)^a;;^ 



{26i — 2bcosuJoT + 2cujosinuJoT) —a]. (20) 
Indeed P and Q can be rewritten as: 
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:[{bk2 + 3A6 + 3c7cjq) coswqt- 



(ca;o/i;2 + 3cAa;o + 375a;o) sinwor — Sjauj^ — {k2 + 3A)(5i], (21) 



a;g[97^a;5 + (^2 + 3A)^J 
8u}o{c6i + ab) sinwoT— 8{b6i + couq) cosojqt], (22) 

with 

K = Aoj^a^ + Auj^c^ + A6j + Ab'^ - . (23) 



Now, solving the equation ( 18 1, we obtain 



Ao = ^P± VP^-Q (24) 
with the conditions P > 0,Q > 0, P'^ > Q or Q < 0. 

We denoted by Aqc and Aqu the amplitude of the controlled system (i.e. when 
6/0 and c / 0) and the amplitude of the uncontrolled system (i.e. when 6 = 
and c = 0) respectively. In the boundary of parameter space where the control of 



the amplitude is efficient, Aqu = Aqc [16 J . Using these conditions, from eq.(24i, 
we find the boundary separating the domain where the control is efficient from the 
domain where the control is inefficient which is given by 



flc = ^p^^^[{^u;la^ + ^Sl){P,-Puf-Qi{P^ + PcPu) 

+ u;o'[97H' + (fe2 + 3A)2]J' ^^'^ 
where Qi is defined by 

Qi = Ab'^ + 4u>qC — 8oJobcsm2u}oT + 8uJo{c6i + ab) sin uJoT 

—8{b6i + cawg) cos wqt, (26) 
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where Pc and P„ represent the function corresponding to the controlled and 
uncontrolled system. 

To validate our technical control in reducing the amplitude vibration, we have 



simulated numerically the set of eq.(25 1. We have plotted in Fig.l for two condi- 
tions c> b and b > c, the domain in the space parameters(T, Oc) where the control 
is efficient in reducing the amplitude of the oscillations, (a) represents the case 
where k2 = and (b) corresponds to A;2 7^ 0. We remark that it is clear that the 
critical value of amplitude excitation a is lower in the case where the control gain 
parameters verify the inequality c > b. We therefore conclude that this domain 
is reduces with k2. The Fig. 2 represents the frequency- and force-response curves 
resulting from eq.(|24]) and the results obtained in Fig. 1 are clearly verified. We 



see that the steady-state response of the system is lower when the parameter a is 
below the boundary domain obtained in Fig. 1(b). Now we have plotted in Fig.3 
the surface amplitude for two conditions for b, where c is fixed, (a) corresponds to 
b = 0.3, c = 0.5 while for (b), b = 0.7, c = 0.5. We noticed from Fig.3(a) that 
the maximum value of the amplitude of steady-state solutions are lower compared 
to the one in Fig.3(b). This differents simulations prove that our technical are con- 
firmed. From Fig.3, one can see that there are range of time delay r and detuning 
parameters 6i where the steady-state solution is zero. This can be a consequence 
of noise effect. 




Figure 1 : Boundary criterion for the effectiveness of control of the oscillation am- 
plitude in the space(ac, r) for (a) k2 = and (b) ^2 = 0.1; blue curve correspond- 
ing to 6 = 0.7 and green curve corresponding to 6 = 0.3, the other parameters are 
7 = 0.002; Q = -0.002; X = -l;5i= 0.01; c = 0.5. 
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Figure 2: The effect of reducing the steady-state response with control parameter 
a obtained in Fig.l (b). (a) Frequency -response and (b) force-response. 




Figure 3: Time-delayed position gain h effect on response-surface curves given by 
Eq. (22): left curve corresponding to 6 = 0.3 and right curve corresponding to 
h = 0.7. 
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3 Effects of the parameters in presence 



In this section, we study separately the effects of gain parameters, parametric ex- 
citation, and combined gain parameters and parametric excitation in the response 



from eq.(24 1 of our model amplitude. 



3.1 Effects of parametric excitation alone with null value of r 



The effects of parametric excitation are found from eq.(24 1 with times delay equal 
zero. We have plotted in Fig.4, the frequency- and forces-responses of the non 
trivial solution for 6 7^ and c 7^ 0. For each values chosen for amplitude or 
frequency of parametric excitation, we plotted two curves. In Fig.4(a) , we see the 
limits of curves correspond to stable and unstable for the nontrivial solution. For 
each value of a, the two curves are the same from certains values of 61. We noticed 
also that for fixed parametric amplitude, two symetric branches appear around the 
point near the origin. The reduction effect due to the parameter a is quite important 
and is found to be visible only in a small region of detuning parameter surrounding 
this point. The Fig.4(b) represents the force-response for two differents values of 
5i. It is clear that the Hopf bifurcation with multi-solution occurs as the detun- 
ing parameter 61 increases. We have also found that the saddle-node bifurcation 
disappears for lower values of 61. 




Ao, 



' 1 


' 1 ' 1 ' 1 
















s, , 


1 . 1 1 , 1 





4-2 O 2-4 6 S 10 



(a) 



(b) 



Figure 4: (a) The effect of the parametric amplitude in the frequency-response of 
the system in the absence of the time-delay, (b) Force-response for two detuning 
parameter values. Here h = 0.3, 7 = 0.002, a = -0.002, A = -1, k2 = 0.1, and 
c = 0.5. 
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3.2 Effects of control gain parameters for a equal to zero 

In this part we have found the control gain parameters effects when the parametric 
excitation a equal to zero. This Fig.5 and Fig.6 illustrated the frequency-response 
of the system for different situations for fixed parameters, (a) and (b) represent 
for each figure the frequency-response of system respectively for a = 0, c = 
and a = 0, c = 0. In Fig.5 r = 0.008,7 = 0.8, a = -0.8, A;2 = 0.1 and Fig.6 
r = 0.008, 7 = 0.002, a = -0.002, k2 = 0.8. We obtain in each case of this 
figure that unstable and stable solutions curves are the same. We have also found 
that the two control gain parameters increases the region of the detuning parameter 
which is more visible when the time-delayed position is considered alone than 
the case where the time-delayed velocity is considered alone. Fig.5 compared to 
Fig.6, prove that the self-excitation parameters a, 7 and ^2 affected also the peak 
of amplitude of vibration which is greater in the case of time-delayed position. 




-2 O 2 4 6 8 10 -2 O 2 4 6 8 10 



5 1 S, (b) 

Figure 5 : (a) The effect of the position time-delay alone in the frequency -response 
of the system without parametric excitation and (b) the effect of the velocity time- 
delay alone in the frequency-response of the system without parametric excitation 
for T = 0.008,7 = 0.8, a = -0.8 and k2 = 0.1. 

3.3 Effects of combined parametric excitation and control gain pa- 
rameters 

In this part, we have plotted in space the (^O; ^i) and {Aq, a) corresponding re- 
spectively to frequency-response and force-response for non-zero time delay , 7 = 
0.8, a = —0.8, k2 = 0.8. This Fig.7 are used to analysis the effects of combined 
parametric excitation and control gain parameters. We have developped the simi- 
lar comments obtained in Fig.4 but the domain of differents bifurcations is not the 
same. 
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Figure 6: (a) The effect of the position time-delay alone in the frequency -response 
of the system without parametric excitation and (b) the effect of the velocity time- 
delay alone in the frequency-response of the system without parametric excitation 
for T = 0.008, 7 = 0.002, a = -0.002 and /c2 = 0.8. 




Figure 7: Effect of parametric perturbation added to the time-delayed term in the 
response of the system, (a) Frequency-response amplitude and (b) Force-response. 
The values of the parameters are r = 0.008, 7 = 0.8, a = —0.8, k2 = 0.8, h = 
0.3, A = -1, andc = 0.5. 
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4 Stability analysis of solutions 



In this case, consider the equilibrium points defined by eq.(24i. To study the be- 
havior of the steady state, we resort to the linearized stability principle by using the 
Routh-Hurwitz criterion lITSll . We apply a little perturbation (with u « Aq and 

V « (t>o ): 



A = Aq + u 

(j)= (j)Q+V 



(27) 



We obtain eigenvalues equation of the Jacobian matrix defined at {Aq, (/)q) of the 



corresponding system equations ( 16 1 and( 17 1 



f]^ + Tr] + D = 0, 



(28) 



where T is the opposite of the trace of the Jacobian matrix and D the determinant 
of the Jacobian matrix. T and D are defined by: 

1 „ 

T = -[-Aujoa + 37a;o^o(l + -4o) - AbsinujoT + Acloq cos wqt], (29) 



D = AoDo = Ao{k2A^ + KiAl + kq), 



(30) 
(31) 



Ki = X + 4wo(367Wo + ca;Q(A;2 + 3A)) sincjoT 
—ijjQ{9c-fU}Q + Ab{k2 + 3A)) coscjJqt, 



(32) 



and 



kq = K — Sbcujo sin 2ujot — 8{bSi + cawg) cos ujqt 
+8(6aa;o + cujo5i) sinwoT 

X = SooQa'y + 97aa;Q + 4:5i{k2 + 3A)a;o — Scjuiq. 



Equation (28 1 has in general two roots: 



-(-r± -4D) 



(33) 



(34) 



(35) 



A positive real root indicates an unstable solution, whereas if the real parts of the 
eigenvalues are all negative then the steady-state solution is stable. When the real 
part of an eigenvalue is zero, a bifurcation occurs. A change from complex roots 
with negative real parts to complex roots with positive real parts would indicate the 
presence of a supercritical or subcritical Hopf bifurcation. The question of which 
possibility actually occurs depends on the nonlinear terms. The Routh-Hurwitz 
criterion implies that the steady-state response is asymptotically stable if and only 
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if T > and D > (i.e.Do > because amplitude Aq are positive) which keep 
the real parts of eigenvalues negative. We obtained 

45sina;oT + 4ca;o coscjot > 4a;oa — 3'JUJqAo{1 + Ao) (36) 

= AujQa'^ + 461 + + 4woC^ - 8c6a;o sin 2ujot- 
8{b5i + cawg) cosa;or + 8(6aa;o + cloq6i) sinu;o''"+ 
K2AI + KiAl (37) 

From these conditions, when parametric excitation and time-delayed position 
(resp. time-delayed velocity) are combined, the stability domain in space {61, a) 
of the non-trivial solutions Aq is obtained in Fig. 8, where (a) correspond to the 
case b < c and (b) correspond to the case b > c. We noticed clearly that for each 
case, the plane is divided in three regions where one of them is unstable and two 
are stable. Other remark is that the stable domain increases when the time-delay 
control gain position is greater than the time-delay control gain velocity. These 
two main conclusions confirm our results obtain in the linear control part and in 
the effect of differents part. 




Figure 8: Stability of non-trivial solution in the space {61, a) with a = —0.1, 7 = 
0.1, c = 0.5, k2 = 0.1, and r = vr. (a) 6 = 0.3 and (b) b = 0.7 
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5 Conclusion 



In this paper we have investgated the interaction of parametric excitation with 
time-delayed position and velocity on the parametric resonance. By the averag- 
ing method to first order, we found the equilibrium point of the system giving the 
amplitude and phase of vibration of the system starting. This study shows the area 
in which the amplitude of the vibration of the oscillator is reduced i.e. the area 
where vibration control is effective. We obtain also for this oscillator the Hopf 
bifurcation or saddle-node bifurcation for certains values of parametric parameters 
and time-delay and we have studied the influence of of parameter k2 which is one 
parameter which modify the ordinary Rayleigh-Duffing oscillator. We applied the 
Routh-Hurwitz criterion for the stability study of steady-state response and we ob- 
tain the stability domain of the parametric oscillator modified Rayleigh-Duffing. 
We have shown clearly that the amplitude of the vibration at the primary resonance 
can therefore be controlled by the active control. Finally, the effects of various 
parameters on the control of parametric resonance have been studied. 

Acknowlegments 

The authors thank IMSP-UAC for financial support. 

References 

[1] Nayfeh, A. H., and Mook, D. T., 1979, Nonlinear Oscillations, Wiley, New 
York. 

[2] Hayashi. C, 1964, Nonlinear Oscillations in Physical Systems, McGraw-Hill, 
New York. 

[3] Chedjou,J. C, Fostin, H. B., and Woafo, R, 1997, "Behavior of the van der 
Pol Oscillator with Two External Periodic forces ",Phys. Scr., 55,pp,390-393. 
DEA de Physique des Uquides. Paris VI- Ecole Polytechnique. 

[4] Carrol, T. L., 1995, "Communicating With Use of filtered. Synchronized, 
Chaotic Signals." IEEE Trans, Circuits syst.,I: Fundam. Theory Appl.,42,pp. 
105-110. 

[5] T. Aida, K. Kawazoe, S. Toda, Vibration control of plates by plate-type dy- 
namic vibration absorbers, J. Vib. Acoust. 117 (1995) 332-338. 

[6] Y. Okada, K. Matsuda, H. Hashitani, Self-sensing active vibration control 
using the moving-coil-type actuator, J. Vib. Acoust. 1 17 (1995) 41 1^15. 

[7] K. Hackl, C.Y Yang, A.H.-D.Cheng, Stability, bifurcation and chaos of non- 
linear structures with control — 1. Autonomous case. Int. J. Non-Linear Mech. 
28 (1993) 441^54. 



13 



[8] A.H.D. Cheng, C.Y. Yang, K. Hackl, M.J. Chajes, Stability, bifurcation and 
chaos of non-linear structures with control — II. Nan autonomous case, Int. J. 
Non-Linear Mech. 28 (1993) 549-565. 

[9] Chyuan-Yow Tseng, Pi-Cheng Tung, Dynamics of a exible beam with active 
non-linear magnetic force, J. Vib. Acoust. 120 (1998) 39-46. 

[10] Aubin, K., Zalalutdinov, M., Alan, T, Reichenbach, R. B., Rand, R. H., 
Zehnder, A., Parpia, J. and Craighead, H. G., Limit Cycle Oscillations in 
CW Laser-Driven NEMS, Journal of Micro-electrical mechanical System 
13:1018-1026,2004. 

[11] Zalalutdinov, M., Olkhovets, A., Zehnder, A., lUc, B., Czaplewski, D. and 
Craighead, H. G., Optically pumped parametric amplification for micro- 
mechanical systems. Applied Physics Letters 78:3142-3144, 2001. 

[12] Rand, R. H., Ramani, D. V, Keith, W. L. and Cipolla, K. M., The quadrati- 

cally damped Mathieu equation and its application to submarine dynamics. 
Control of Vibration and Noise: New Millennium 61:39-50, 2000. 

[13] Wirlcus, S., Rand, R. H. and Ruina, A., How to pump a swing. The College 
Mathematics Journal 29:266-275, 1998. 

[14] Zhehe Y,Deqing M., Ziehen C, Chatter suppression by parametric exci- 
tation: Model and experiments. Communications in Nonhnear Science and 
Numerical Simulation 330:2995-3005, 2011. 

[15] Xin-ye L., Yu-shu C, Zhi-qiang W. and Tao S., Response of parametrically 
excited Duffing-van der Pol oscillator with delayed feedback , Applied Math- 
ematics and Me- chanics 27:1585-1595, 2006. 

[16] M. Siewe Siewe , C. Tchawoua , S. Rajasekar, Parametric Resonance in the 
Rayleigh-Duffing Oscillator with Time-Delayed Feedback. 

[17] Darya V. Verveyko and Andrey Yu. Verisokin, Application of He's method to 
the modified Rayleigh equation,DiscTete and Continuous Dynamical Systems, 
Supplement 2011, pp. 1423-1431. 

[18] Guckenheimer J., and Holmes P., Nonlinear oscillations, dynamical systems, 
and bifur- cations of vector fields. New York: Springer- Verlag; 1983. 

[19] J.C.Chedjou, L.K.Kana,I.Moussa,K.Kyamakya,A.Laurent, Dynamics of a 
quasiperiodically forced Rayleigh oscillator, TrsnsacMons of the ASME, 600/ 
vol.128, Sptember 2006. 

[20] B.R. Nana Nbendjo, Y Salissou, P. Woafo, Active control with delay of 
catastrophic motion and horseshoes chaos in a single well Duffing oscilla- 
tor. Chaos Solitons and Fractals 23 (2005) 809-816 . 



14 



[21] R. Tchoukuegno, B.R. Nana Nbendjo, R Woafo, Linear feedback and para- 
metric controls of vibrations and chaotic escape in a potential, Interna- 
tional Journal of Non-Linear Mechanics 38 (2003) 531-541. 

[22] Nayfeh A.H., (1981), Introduction to Perturbation Technique, J.Wiley, New 
York. 



15 



